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H ■ In this paper we describe formulas, derived from conjectures of Conrey, Farmer, 

Keating, Rubinstein, and Snaith [CFKRS] . for the moments of quadratic Dirichlet L- 
functions at the central point, and the moments of L-functions associated to quadratic 
twists of an elliptic curve. 

We are motivated to study moments in these two families of L-functions because 
of their apparent connection to the moments of characteristic polynomials of unitary 
symplectic and orthogonal matrices. 

Montgomery was the first to discover a link between an L-function and characteristic 
polynomials of unitary matrices [Moj . He computed, with restrictions on the allowed 
test functions, the limiting pair correlation of the zeros of the Riemann zeta function, 
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and found that it coincides with the average pair correlation of the eigenvalues of large 
random (according to Haar measure) unitary matrices that had been computed earlier by 
Dyson |Dy| . Odlyzko later confirmed this agreement numerically, without restriction [0] . 
Rudnick and Sarnak generalized Montgomery's result to higher correlations and to any 
primitive L- function |RSj . 

Katz and Sarnak then made precise connections between various families of L- functions 
and matrices from specific classical compact groups, based on results linking the density 
of zeros L-functions and analogous zeta functions over function fields, to the eigenvalue 
densities of random matrices in the classical compact groups |KSj |KS2] . For instance, 
their work showed a statistical connection between the zeros of quadratic Dirichlet L- 
functions and eigenvalues of unitary symplectic matrices, and between the zeros of L- 
functions of quadratic twists of an elliptic curve and eigenvalues of orthogonal matrices. 
The papers [Rj and |R2] . provided further theoretical and numerical support for the 
relevance of these matrix groups to our two families of L-functions. 

Subsequently, Keating and Snaith were able to predict the leading term in the asymp- 
totics for the moments of the Riemann zeta function on the critical line by carrying out 
an analogous computation for random unitary matrices [KeSj . In a companion pa- 
per [KeS2J . they also conjectured the leading term in the asymptotics for the moments 
in our two families of L-functions by computing the moments of the characteristic poly- 
nomials of random unitary symplectic and even orthogonal matrices. See also the paper 
of Conrey and Farmer |CF| which contains some arithmetic information needed for the 
Keating and Snaith approach to moments. 

The method of Keating and Snaith for predicting moments of L-functions relies on 
computations in random matrix theory, for example it uses Weyl's integration formula 
and the Selberg integral, and some guesswork to make the heuristic leap to number 
theoretic moments. It also has the drawback of requiring, as input, the relevant classical 
compact group as predicted by Katz and Sarnak. 

On the other hand, the approach, referred to above, of Conrey, Farmer, Keating, 
Rubinstein, and Snaith does not rely on random matrix theory to derive, heuristically, 
the moments of various families of L-functions. Their method is based strictly on num- 
ber theoretic techniques involving the approximate functional equation, the traditional 
equation that is used to study moments of L-functions [T] [j]. While random matrix 
theory is not needed in their approach, the formulas that their heuristic approach yields 
for L-functions have provable analogues in random matrix theory. CFKRS were also 
able to make progress by introducing 'shifts' into the moments, a strategy that was 
inspired by Motohashi's evaluation of the fourth moment of the zeta function |Motj 
and also by an analogous problem in random matrix theory. Their method, therefore, 
produces an answer that can be compared against various moment computations in ran- 
dom matrix theory, and, instead of using the predictions of Katz and Sarnak, it provides 
evidence for them. Furthermore, the conjectured formulas of CFKRS go well beyond 
the leading asymptotic of Keating and Snaith, providing, implicitly, a full asymptotic 
expansion for a variety of L- function moment problems. Because their conjectured for- 
mulas provide a full asymptotic expansion for moments, one can test them numerically 
by comparing the predicted moments against those computed from L-function data. See 
for instance |CFKRSj |CFKRS2j [SR] [RY] . 
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Our goal is to turn the implicit formulas of CFKRS into asymptotic expansions with 
explicitly given coefficients. We elaborate on the CFKRS formulas, for the family of qua- 
dratic Dirichlet L-functions, in Section fl.il and for quadratic elliptic curve L-functions 
in Section [5l 

Besides the approaches of Keating and Snaith and of CFKRS, two additional methods 
have yielded interesting results for the moments of L-functions. 

Gonek, Hughes, and Keating [GHKj . and Bui and Keating [BK] use the explicit 
formula for an L- function to realize the L- function as a hybrid between partial Hadamard 
and Euler products. They assume statistical independence between these two products 
and study the moments of the partial Euler product using number theoretic heuristics. 
The moments of the partial Hadamard product are studied by modeling the zeros of 
the Hadamard product based on the predicted classical compact group. Their approach 
therefore suffers the same disadvantage of the Keating and Snaith method of requiring 
the predictions of Katz and Sarnak as input. The main advantage of their method over 
the Keating and Snaith method is that it explains, rather than guesses, the appearance 
of an 'arithmetic factor' in moment formulas for L-functions. Another disadvantage is 
that it only seems to correctly predict the leading asymptotic for the L-function moments 
that they consider, and thus only agrees with the CFKRS prediction to leading order. 
Presumably this is because their assumptions are too strong, for example the statistical 
independence between the partial Hadamard and Euler products, and their use of matrix 
eigenvalues to model the partial Hadamard product. 

Another method for studying moments of L-functions has been developed by Diaconu, 
Goldfeld, and Hoffstein |GHj [DGH| and uses the theory of multiple Dirichlet series. It 
has the advantage of proving asymptotic formulas for some L-function moments, for 
example the first three moments of quadratic Dirichlet L-functions at the central point. 
However, it has the disadvantage of involving an elaborate sieving process (in the case 
of quadratic characters), that makes it unwieldy for producing explicit formulas for the 
asymptotic expansion. Interestingly, their method predicts the existence of additional 
lower order terms of smaller magnitude that go beyond those of the asymptotic ex- 
pansion of CFKRS. See the paper of DGH as well as that of Zhang [Z], and Alderson 
and Rubinstein |AR| for discussions and computations regarding these additional lower 
terms. 

1.1. The CFKRS conjecture for L(l/2,Xd). We begin by describing the CFKRS 
conjecture for quadratic Dirichlet L-functions. Let L* be a squarefree integer, L* 7^ 0, 1, 
and let K = Q(vi^) be the corresponding quadratic field. The fundamental discriminant 
doi K equals D if D = 1 mod 4, and 4L> if L> = 2, 3 mod 4. Let Xdin) be the Kronecker 
symbol f ^j, and L{s, Xd) the quadratic Dirichlet L-function given by the Dirichlet series 

00 / \ 

(1.1) L(5,Xd) = E^' 3f?(s)>0, 

satisfying the functional equation 

(1.2) L{s,Xd) = Idl'^-' X{s,a)L{l - s,Xd) , 
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where 

1— s+a 



('•') X(.,a)=vr ^^^(ip^, «=|^ .f,^0. 

Let S{X) denote the set of fundamental discriminants with \d\ < X. The Gamma 
factor in functional equation for L{s, Xd) depends on whether d < or d > 0. Thus, 
define further 

S+{X) = {d€ S{X):d>0} 
(1.4) 5_(X) = {deS{X):d<0}, 

to be, respectively, the sets of positive and negative fundamental discriminants with 
\d\<X. 

CFKRS conjectured [CFKRS] the asymptotic expansion: 

(1-5) Y. L{l/2,Xd)''--^XQ±{k,\ogX), 

— TT 

des±{X) 

where Q+{k,x) and Q-{k,x) are polynomials of degree k{k + l)/2 in x that we will 
describe below. The fraction S/vr^ accounts for the density of fundamental discriminants 
amongst all the integers. 

The polynomial Q± (k, log X) is expressed in terms of a more fundamental polynomial 
Q±{k,x) of the same degree that captures the moments locally: 

X 



(1.6) Q±{k,logX) = ^J Q±{k,logt)dt 



One of the main achievements of CFKRS was to give a general recipe/heuristic for 

producing formulas for moments of various families of L-functions. Their formula (see 

Conjecture 1.5.3 in [CFKRS ] ) for the polynomial Q±{k,x) is given implicitly in terms 

of a A;-fold multivariate residue: 

(1.7) 

^ /, ^ (-l)'=(^-i)/22'= 1 / rG±{zi,...,Zk)A{zl...,zlf ^y-fe ,, , 

Q.iK^) = '-^^^ (2^/-/ ntJ.f-' '^ ^'^-''^^^-^^^^ 

where A(u;i, . . . ,Wk) is the Vandermonde determinant 

(1.8) A{wi,...,Wk) = det{wl~ )kxk= Yl i'^-i - ^i) 



l<i<j<k 



and 



k 

(1.9) G±{zi,...,Zk) = Ak{zi,...,Zk)f[X{l + z,,ar'/^ [] C(l + ^i + 



j=l - l<i<j<k 



Here, a = for G+ and a = 1 for G_, X{s, a) is given in (jl.3|) . and A^ equals the Euler 
product, absolutely convergent in a neighbourhood of (zi, . . . , Zk) = (0, . . . , 0), defined 
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by 

(1.10) Ak{z,,...,zk)=i[ n (i 



jjl-j-Zi-f-Zj 

P i<i<j<k 




One advantage of equation (jl.7p is that it allows one to easily see that Q±{k,x) is a 
polynomial of degree k{k+l)/2 in x. That is because the denominator of the multivariate 
residue picks up terms in the numerator involving n^=i ^j ~ ! which is of degree 2k{k — 
1). Now, the factor A^zf, . . . , z^)^ is a homogeneous polynomial, also of degree 2k{k — l). 
However, the factor G±{zi, . . . , Zk) cancels k{k + l)/2 of the factors of the Vandermonde, 
because each ^(1 + Zi + Zj) has a Laurent expansion that begins l/{zi + Zj) coming from 
the pole at s = 1 of C(s). Therefore, in considering the multivariate Taylor expansion 
of the numerator about zi = . . . = Zk = 0, we only need to take terms in the series 

(1-11) '"p( 2^^/^- ) =E^(-i + --- + -^r 

up to n = k{k + l)/2. Hence, in the x aspect, the /c-fold residue only involves terms up 

to x^^^+^y\ 

Equation (|1.7|) has the disadvantage of expressing Q±{k, x) implicitly. Let us therefore 
write 
(1.12) Q±{k, x) = c±(0, A:)x'=('=+i)/2 ^ ^^(^^ f^^^k(k+i)/2-i ^ _ _ _ ^ ^^^^f^ ^ ^y^, k). 

Our main result, described in the following theorem, gives explicit formulae for the 
coefficients c±{r,k). We first define 

fc(fc+i) / , -, \-fc / -, \ -fc 



(1.13) a.:=^.(0.....0)=n ''7il' T' ^ "" ^ "^ +' 

Theorem 1.1. In /il.l2\) . the leading coefficient c±{0,k) of Q+{k,x) or Q-{k,x) are 
both equal to 

and, for given r > 1, we have 

(1.15) c±{r, k) = c(0, k) J2 bfik)Nxik), 

|A|=r 

where Nx{k) is a polynomial in k of degree at most 2\X\, defined in (|2.52|) . a^ is defined 
in ()1.13p . and the b^ {k) 's are the Taylor coefficients of a holomorphic function, defined 
in (j2.4p and (j2.5p . The sum is over all partitions \X\ = r, with ^ Aj = r and Ai > A2 > 
... >0. 
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We remark that formula (|1.14p for the leading term agrees with the prediction of 
Keating and Snaith. See (34), (45), and (47) of Keating and Snaith [KS2j (replacing 
logD by X in their equation (45)). Their derivation is heuristic and based on the 
Selberg integral. Compare also to the leading term of equation (1.5.17) of [CFKRSj . 
with N = x/2 in that equation. To verify the agreement between these, one can check, 
inductively, that: 

1 ''~^ (2j)! 1 '^ 1 '^ .-I 



/i (^ + JV- 2M^'+i)/2 11 (2j - 1)!! 11 (2j)! 



Note that (jl.lSp is analogous to formula (1.16) of jCFKRS j which provides a formula 
for the coefficients of the moment polynomials of the Riemann zeta function. See also 
Dehaye's paper [D], also for the Riemann zeta function, where he gives a combinatorial 
formula for the analogue of our polynomial Nx{k). 

We work out examples, for r = 1 and r = 2. Table [1] provides Nn\{k) = k(k + 1), 
N^i^i){k) = \k{k - l){k + l){k + 2), and iV(2)(A;) = 0. Thus, 

c±(l,A;) = c(0,A:)6j)(A:)iV(i)(fc) 



fc-i 



(1-17) =|n(if7yTMfc + i)^(i)(^)- 



and 



c±(2,fe) = c(0,A;) (6±^^(A:)iV(i,i)(A;) + 6f2)(A:)iV(2)(A:)j 

(1-18) = I n (if7)T X ^(^ - 1)(^ + 1)(^ + 2) ^fi,i)(^)- 

Let 



-. oo 

(1-19) c(i + ^) = - + E(-ir7"^' 



■s ' — ' n\ 



n=0 

be the Laurent expansion about of C(l + s) (70 is Euler's constant), and define 

(_l)i(pl/2 _ lyj-k + (pl/2 ^ ^yj-k 



(1.20) f,{p) :-- 



(pi/2 _ i)-fc + (pi/2 + lyk + 2p-i-fc/2 ' 



Formulas for the coefficients b,-^Ak), b,^-^Ak) can be derived using the method de- 
scribed in Section [3l and are given by 

(1.21) 6±^(A;) = -llogvr + i^(l/4 + a/2) + (fc + l)7o + E(^^ + /i(p))log^' 
where a = for 6^w and a = 1 for bV^., and 

(1.22) 6±^(A;) = hl^{kf-^l-2^,-Y.(^j-^ + h{pf-h{p)yog{pf. 



LOWER ORDER TERMS 



In Section[2]we derive formula ()1.14p for the leading coefficient of Q±{k, x). Our tools 
are then applied, in Section [2.21 to the general term c±{r, k), where we obtain a formula 
for Nx{k) expressed as a sum of determinants of the form: 

2k — i — Afc_j-|-i 



(1.23) Dx{k) = det^^ 

/ / l<i,j<k 

where A = (Ai, . . . , A^) is a partition with length 1{X) < k (see Section 11.21 for defini- 
tions) . 

In Section |4] we derive some interesting formulas for these determinants. To describe 
our formulas, let y = {yi, . . . ,ym)- We define the coefficient operator [yf^] on the set 
of formal multivariate Taylor or Laurent series in y, which picks the coefficient of the 
monomial y" in the series. More precisely, if 

(1-24) /(yi,...,y„)= Y, ar„...,rrny?---y"r;r, 

ri,...,rm,eZ 

define 

(1-25) [yT---ynr]f = au,,...,u^. 

We prove the following Theorem. 

Theorem 1.2. Let X be a partition and fi be the conjugate partition. Let m = 1{X), and 
n = l[^) = Ai. For k > max(/(A), Ai), we have 

(1.26) Dxik) = 2(2)-!^! 

[yt^"'-'...yt]{ n (y-yi)(i-y-y.) 11(1-2/')-'"" 



X 



l<i<j<m 1=1 



and also 

(1.27) Dx{k) = 2(2)-!^! 



X [z('^+""-^ . . . z^"l I TT r^. _ r.vi ^ ^. ^ r.^TTn -L9^,Vi ^ ^.^fc-** 



-1 






Corollary 1.3. Let X be a partition, with /(A) = m. There is a polynomial P\{h), 
integer valued at integers, of degree \X\ such that for k > max(/(A), Ai), 

(1.28) Dx{k) = 2(2)-!^! X Px{k). 
The leading coefficient of Px{k) is 

(1.29) ^l,.<.<n.(A.-A,-^ + ,^) ^ 

]\l<i<m{h + m - l)\ 

where x'^(l) is the degree of the irreducible representation of the symmetric group S\x\ 
indexed by X. In particular, 

(1.30) Z)o(A;) = 2(2). 

where DQ{k) is the determinant associated to the empty partition. 
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Table [2] gives a list of the polynomials Px{k) for partitions up to weight 7. Observe, 
in the table, that Px{k) often has many linear factors. This fact plays a role in our 
formula for Nx{k) so we encode it in the following corollary. 

Corollary 1.4. Let X be a partition. Then Px{k) is divisible by 

(1.31) {k - Xi){k - Ai - 1) ... (A; - /(A) + 1) x {k + Ai)(A; + Ai - 1) . . . (A; + 1{X)), 

where we take the first product to be 1 if Xi > 1{X), and the second product to be 1 if 
Ai </(A). 

Finally, in Section[5]we discuss the application of our techniques to the related problem 
of the moments of the L-functions associated to quadratic twists of an elliptic curve. 

1.2. Symmetric function theory. We collect here some definitions and results from 
the theory of symmetric functions that we use in our paper. The details can be found 
in [Ml Chapter 1]. We have used the notations of [M] . 

A partition A is a sequence of non negative integers (Ai, A2, • • • ) such that 

(1.32) Ai>A2>---, 

and only finitely many AjS are non zero. 

The length of the partition A is defined to be the number of non zero AjS. We denote 
it by 1{X). The weight of a partition A, denoted by |A| is 

(1.33) |A|=EA- 

The diagram of a partition is the set of points 

(1.34) {(^,j)|l<^</(A), 1 < i < AJ. 

The conjugate partition X' of a partition A is the partition whose diagram is 

(1.35) {{hi) I (ii^) is in the diagram of A}. 
Equivalently, the conjugate partition of A is a partition A' = (A'^, A2, ■ ■ ■ ) where 

(1.36) A', = #{A, I A,- > i]. 

The symmetric group Sn acts on the polynomial ring Z[xi, . . . , x„] by permuting the 
independent variables xi, . . . , x„. The ring of symmetric polynomials in n- variables, A„, 
is the set of polynomials in 1j[xi, . . . ,x„] which are invariant under this action of Sn- 
The ring A„ is a graded ring: 

(1.37) An = 0A^, 

fc>0 

where A^ is the set of homogeneous symmetric polynomials of degree k. 
For m > n, there is a ring homomorphism 

(1.38) pm,n ■ 1'[xi,...,Xm] -> Z[xi, . . . , X„], 

where Prn,n{xi) = Xi for i < n, and Pm,n{xi) = for i > n. This restricts to a map 

(1.39) pra,n ■ A^ ^ A^ 
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The maps given by (11.390 define an inverse system. Let 

(1.40) A'= = ^A^, 
and 

(1.41) A = 0A^ 

A:>0 

The ring A is cahed the ring of symmetric functions. This is a graded ring. The definition 
of A gives us maps 

(1.42) /On : A ^ A„. 

In this paper, we shall use four Z-bases, parametrized by partitions, of the ring A: 
the monomial symmetric functions {mx), elementary symmetric functions {e\), complete 
symmetric functions (hx) and the Schur symmetric functions {sx)- In addition, we shall 
be using power symmetric functions (px). The power symmetric functions form a Q 
basis of A (g)^ Q. We shall use the same symbols to denote their image under pn in A„. 

Given a = (ai, . . . ,an), we write x" to denote x"^ • • • x"". Let A be a partition of 
length less than or equal to n. We define the monomial symmetric function mx by its 
image under pn for every n. If n > /(A), then 

(1.43) mA(xi,... ,x„) =^x", 

a 

where the a ranges over distinct permutations of (Ai,...,A„). If /(A) > n, then 
?n,A(xi, . . . , Xn) = 0. For the only partition of 0, the empty partition, we define m^ = 1. 
Let r > be an integer. The elementary symmetric function e^ G A is given by 

(1.44) er= ^ Xi^ ...Xi^ = m(i_...^i), 

and eo = 1. For a partition A, we define 

(1.45) ex = ex^ex2 ■ ■ ■ ■ 
The generating function for Cr is 

(1.46) E{t) = Y,erf = W{l + Xit). 

Let r > be an integer. The complete symmetric function h^ is defined to be 

(1.47) hr=J2 "^A. 

|A|=r 

Given a partition A, we define 

(1.48) hx = hx.hx, . . . . 
The generating function for hr is 

(1.49) H{t) = Y,hrf = Y[{1-Xity\ 

r>0 i>l 

Equations (jl.46p and (|1.49p give us the identity, 

(1.50) H{t)E{-t) = 1. 
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For r > 1, the power symmetric function pr is defined as 

(1.51) pr = J2^l = '^ir)- 

i>l 

For a partition A, we define 

(1.52) PA = PA1PA2 • • • • 

Let (ai, . . . , a„) G N". We define a^ G Z[xi, . . . , x„] by 

(1.53) aaixi,... ,Xn) = det{x"')i<ij<n. 

Clearly a^ is skew-symmetric; that is, for w G Sn, w{aa) = sgn{w)aa, where sgn(^i;) is 
the sign of permutation w. Let Sn be the partition 

(1.54) 5„ = (n-l,n-2,...,l,0). 

For a partition A of length less than or equal to n, we append O's as necessary to A to 
create an n-tuple, and define 

(1.55) Sx[Xl,...,Xn) = -, r-. 

0'S„[Xl, • • • , Xn) 

This is a polynomial. Since sx{xi, . . . ,Xn) is a ratio of skew-symmetric polynomials, 
it is a symmetric polynomial. These symmetric polynomials are called Schur symmet- 
ric polynomials. For m > n, pm,n{s\{xi, . . . ,Xm)) = sx{xi, . . . ,Xn), hence they are 
represented by a function s\ £ A. 

Using the definitions of ax and sx, it is easy to check that 

(1.56) as„{xi,...,Xn) = Yi i^i-Xj)^ 

l<i<j<n 

and 

(1.57) Ss„ixi,...,Xn) = Yl i^i+Xj). 

l<i<j<n 

Let A be a partition and A' be the conjugate partition. Then, for n > /(A) [Ml P-41] 
(1-58) sx = det (/iA,-i+i)i<ij<„ , 

and for m > /(A') 



(1-59) sx = det[ey^,^,j^^^^^^^ 

Identity (jl.SSp is called the Jacobi-Trudi identity, and ()1.59p is called the dual Jacobi- 
Trudi identity. Schur symmetric functions satisfy [Ml p. 63] 

(1-60) n (1 - ""iyj^ = E ^A(x)sA(y), 

and 

(1.61) n(i + ^^yi) = E^A(x)sA'(y). 

*J>1 A 

The sum in (ll.60p and ()1.6ip is over all partitions A. Identity ()1.60p is called the Cauchy 
identity, and (|1.6ip is called the dual Cauchy identity. 
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There is a fundamental involution to, a ring automorphism, defined on the ring of 
symmetric functions: 

(1.62) u{er) = hr. 
Using (jl.SOp . we can prove that 

(1.63) U}{hr) = Cr- 

We also have 

(1.64) uj{sx) = sy, and a;(p„) = (-1)"" V- 

2. Terms of the asymptotic expansion 

We begin by rewriting the integrand on the right hand side of (jl.7p as a ratio of a 
holomorphic function and a monomial. The function G{zi, . . . ,Zk) in (|1.9p has a pole 
in each Zj at (0, . . . , 0) coming from the product of the zeta functions. These poles are 
eliminated by a portion of the Vandermonde determinants. Note that 

(2.1) A{zl...,zlf=[ n (,^ + ,^.)»^(^i'---'^^)^(^i'---'^i) 



l<i<j<k 



^'uU 



Specifically each factor (zj + Zj) occurring here cancels a pole coming from (^{1 + Zi + Zj] 
We obtain ()2.ip by observing 

A{zl ...,zl) = Hiz] - zf) = A(zi, ...,Zk)X{{z^ + Zj) 

i<j j>i 

(2.2) =A{z,,...,Zk) - , . 

^ i.lj=i^j 

Substituting (|2.ip into (jl.7p . we have 

(_l\k(k~l)/2 Iff 

(2.3) Qi(/c,x) = ^^^^ .^^j...JA,{zu...,z,) 

k 
[] X{\ + Zj,a)-^2 11 (zi + Zj)Cil + Zi + Zj) 

j=l l<i<j<k 

Aizi,...,Zk)A{zl...,zl) fxA \ 



nt=i#-^nt=i-. 



dzi . . . dzk ■ 



Now the integrand is written as a ratio of a function which is holomorphic in a neigh- 
bourhood of (0, ... , 0) and a monomial. 

Recall that a^ = ^^(0, ...,0). Let z = (zi,...,Zfc) and m\{z) be the monomial 
symmetric polynomial defined in (ll.43p . Let 

oo 

(2.4) E E &A (^)"^a(^) 

i=0 \\\=i 
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be the power series expansion of 

k 



(2.5) —Ak{zi,...,Zk)X{X{\ + Zj,a) 2 \{ {zi + Zj)(:{l + z^ + 

"^ i=l l<i<j<k 



Here, the coefficients h~^ are associated to the a = 1 case, and 6^ with a = —1. 



In (|2.4p . the sum is over all partitions Ai + . . . + A^ = i, with Ai > A2 > • • • Afc > 0. 
We divide the expression by a^ to ensure that the constant term in the power series is 
1. We shall calculate the Taylor series of (j2.5p by calculating the Taylor series of its 
logarithm. This calculation is simpler if the constant term is 1, i.e. 6q (/c) = 1 in ([27 
So (12.31) becomes 



(_'\\k{k-l)/2 a °° r r 

(2.6) Q^{k,x)=^-^^ 7^EE&A(^)f ••/-a(.i,...,..) 

A{zi,...,Zk)A{zl...,zl) fxJ^ \ 

^k_^ ^2k «^p ( 2 e^ ^' ) ■ 

Only finitely many integrals in the sum (|2.6p are nonzero. Each of the integrals in (j2.6p 
picks up the coefficient oi Zi ~ . . . Zj^ ^ in the Taylor expansion of the numerator of the 
corresponding integrand. If degmA(-2i; ■ ■ ■ , Zk) + deg A(zi, . . . ,Zk) + deg A(z^, . . . , z^) > 
deg{zi ~ . . . z^ ~ ), that is |A| > k{k + l)/2, then in the Taylor expansion of the nu- 
merator of ([276]) the coefficient of zj''^'^ . . . z^~'^ is 0. 

Given k, and a A in the sum (j2.6p . the coefficient of the monomial z-^ ~ . . . z^ ~ in 

the Taylor expansion of the numerator of the integrand is a constant, depending on A 

, M^±l)_l\l 

and K, times x 2 1 i_ 



2.1. The leading term. In this section, we shall calculate the leading coefficient of 

fc(fc+i) 
Q±{k, x), i.e. the coefficient c±{0, k) oi x 2 . The calculation will also provide insight 

into how to calculate the lower order terms of Q±{k,x). The leading coefficient is the 

same for Q^(k,x) and for Q-{k,x), and is given in the following proposition. 

Proposition 2.1. The leading coefficient c±(0,k) of Q±{k,x) in \1.12\) is 



p.v) sn "'■" 



The leading term in (J1.12p corresponds to the i = Q term of (j2.6p . In this case there 
is only one integral within the inner summation sign, giving 
(2.8) 



k\{2mf J J rii-i-^ 



k- 
1 ^i 3=1 
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Substituting Uj = xzj/2, simplifying, and then relabeling Uj with Zj, we obtain 
(2.9) 

., fc(fc-l) fc(fc + l) ^ , s ^ / 9 On k 

The presence of the Vandermonde determinants prevents us from separating the in- 
tegrals. However, we apply the following trick to move the Vandermonde determinants 
outside the integral. Introduce new variables xi,. . . ,Xk and consider the more general 
integral 

(2.10) 

Thus, the evaluation of c±(0. A;) boils down to determining /(I, . . . , 1). 

Next, we introduce a partial differential operator which will help us move the Van- 
dermonde determinants outside the integral. Note that for a polynomial P{xi, . . . ,Xk) 
in k variables, we have 

(2.11) P (^— , . . . , — j exp ( J2 ^j^j 1 = ^(^1' ■ ■ ■ , zk) exp ij2 ^i^i 

We set 

(2.12) q{zi, ...,Zk):= A(zi, . . . , Zk)A{zf, ..., zl). 
Then (l2lll equals 

^o^Q^ 1 / / (d a \ exp(^|=i^j^ , , 

Pulling the differential operator outside the integral (Leibniz's rule) we conclude that 
(I2J3I1 equals 

d d \ 1 r /exp(^^^iXjZj 



(2-14) gU-.---.17- TTT-^^--'^ '-^-\r' dz,...dz, 



dxi'"" dxk ) (27ri)*= J / J]i=i zf 



The integrand in (|2.14p can be written as a product of integrals in one variable. 



Each integral in the above product can be evaluated by expanding exp(2;_,-2;j) = X^n^ol^j-^j)"/"'' 
.f- is .f 



The coefficient of zf'"^ is xf'^/{2k - 1)!, and thus (l2J5]l equals 



d_ _d_\J^ xf-^ 



We have turned our computation of c± (0, k) into the question of determining the 
result of applying q{-^, • • • , gf-) to ni=i (2k-iY ' ^^^ finding the value of the resulting 
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polynomial at (1, . . . , 1). This calculation is done in Lemma 12.41 The proof of Lemma 

\2H uses Lemmas ESI and EH 

Lemma 12.21 a variant of Lemma 2.1 in [CFKRS2] . gives a formula for applying the 



differential operator A 



92 



Lemma 2.2. 

(2.17) 

Proof. Write 
(2.18) 



A 



2' 



dxf 



A 



92 
9x1 






to a product of functions. 



)k 



(2i-2), 



fcxfc 



5xf ' " ' ' 



^ 

9x2 



a2i- 



.{2i-2) 



fcxfc 



Applying this to Y\i=i f{xi), and noticing that xj only appears in the i-th row of the 
determinant, we can move /(xj) into that row. D 

Lemma 12.31 gives a formula for applying a product of differentials to a determinant of 
functions. 

Lemma 2.3. Let fi{x), . . . , /fc(x) be smooth functions of one variable. Then 



(2.19) 



d" 



dx'l' 



Q-n-k 

dx^^ 



fi{xi 



fl{Xk) 



fk{xi 



fk{Xk) 



p("i) 



("i)/ 



/r(xi) tTi^i) 



f("fc) 



(n-k). 



fr'ixk) fn^k) 



Proof. It is easy to see if we first look at a simple case, say -^r applied to the determinant 
on the left hand side of HUM . 

Lemma 2.4. Let q{zi, . . . , z^) = A(zi, . . . , Zk)A(zf, . . . ,2^), then 



a 



(2.20) 






evaluated at (xi, . . . , Xk) = (1, . . . , 1) is 



(2.21) 



fe(fe-i) 
(-1)^^ X A;! 



'fc-i 

n 



(2j)! 



fe(fc-i) 



Proof. To prove the Lemma, we relate the value of ()2.20p evaluated at (xi, . . . ,Xk) = 
(1, . . . , 1) to a determinant of a matrix whose entries are binomial coefficients. We then 
use an identity for binomial coefficients to rewrite the determinant as a product of two 
determinants, and evaluate each of them separately. 
Applying Lemma [221 we can deduce that 



(2.22) 



d 
dxi 



d \ ^ f d"^ 

A 



dxk 



d^ 



d'^x\ 



2' ■ 



dxi 



7^ n/(^ 



fc/ j=i 
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equals 
(2.23) 



A 



d 



d 



{2(i-l)), 



Xi) 



dxi ' ' dxi^J '' ^' "' kxk 

Expanding the Vandernionde determinant of partial differential operators, we obtain 

/(Xi) /(2)(xi) ••• /(2(^-l))(xi) 



(2.24) 






MGSfe 



dx'l' 



dxl 



f{x2) f'H 



X2 



f{Xk) f'\xk) 



where iJ,i, . . . , fik is the image of the permutation fj, of 1, . . . ,k 
we can see that (|2.24p equals 

/(Mi-i)(a;^) /(mi+i)(^^) .. 



f 



(2(fc-l)) 



{X2] 



(2.25) 



Let fix) 



(2.26) 



^J,eSk 



Applying Lemma [27 



j{|t2-l+2(fc-l))(-^2) 



&'-'Hxk) &^+'Hxk) 



(2fc-l)! 



Expression (j2.25p evaluated at (xi, 



/(A'fc-l+2(fc-l))(xfc) 

,Xfc) = (!,...,!) is 



X] sgn(/i) 



1 


1 


(2fc-Mi)! 


{2k-fii-2y. 


(2fc-M2)! 
1 


(2k-fi2-2y. 
1 



(-M1+2)! 



(-M2+2)! 



(2fc-Atfc)! (2A.-^fc-2)! (-Mfc+2)! 

Rearranging the rows to cancel the effect of /i (this introduces another sgn(/i) in front 
of the determinant) and evaluating at (xi, . . . , x^) = (!,..., 1), we get (|2.26|) equals 



(2.27) 



k\ 



1 



(2fc-l)! (2fc-3)! 
(2Jk-2)! (2fc-4)! 



J_ 1 

fc! (fe-1)! 







We can convert the determinant (j2.27p into a determinant of matrices whose entries are 
binomial coefficients. Multiplying the j*^ column by ,„, ■_^^^| and the i*^ row by {2k — i)\, 
we see that ()2.27p equals 

0!2!---(2fc-2)! 



(2.28) 



A;! 



(2fc-l)!(2A:-2)!---A;! 



f2k-i\ /2fc-i\ 

V / \ 2 ) 

(2k-2\ /2k-2\ 

\ Q I \ 2 ) 



(2) 



\2k-2) 
(2k-2\ 
\2k-2) 

\2k-2) 



The determinant in (|2.28p is 
(2.29) 



2k -i 
2j-2 



kxk 
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In Section Owe study this determinant. From (11.230 and Corollary II. 3^ the determinant 
of this matrix equals (— 2)*^2J. The extra (—lY^) here comes about from the fact that 
the D{)[k) in (jl.23p has its columns reversed from the above determinant. 

D 

Applying Lemma 12.41 to (j2.16p . we find that the leading term is: 



aUx\''T' ( m\■■■{2k- 
k\ \2) \ ■ (2fc-l)!-- 


-2)\\ fc(fc-i) 
.k\) 


(2.30) -''' n\/^^^U'^'+'^^' 




Hence the coefficient of the leading term is 




(2 31) "^ff ^'^'^■ 





This proves Proposition 12. 1[ 

2.2. Further lower order terms. In this section we consider a general term occurring 
in the sum of integrals (12.60 . Let A be a partition. We shall calculate 



(K) Y ■■■ Y inx{^ 

2 .v2^ /^ k 



(2.32) '—L^^ _A_6±(fc)f ..<^m;,(zi,...,zfc 



kl {2Triy 



A(zi,...,Zfc)A(zf,...,z|) Xy;-^ 1 

exp -2^Zj dzi ...dzk- 



2^ nti ^f V 2 f 



i=i 



Modifying the approach of the previous section to incorporate the extra monomial 

m\{zi, . . . , Zfc)) we define 

(2.33) qx{zi,...,Zk) = mx{zi, . . . ,Zk)/^{zi, . . . , Zk)A{zl, . . . ,zl). 

Following the same steps as in the evaluation of the leading term, expression (j2.32p 

becomes 

(2.34) 

(-1) 2 akbf[k) /x\ 2 l-^l / /_d_ _d_\ ^ x- 



k\ \2J \^^\dxi"''dxj^}-^{2k-l)\ 

This section is devoted to calculating ()2.34p . 

Let f{x) = x'^^-^/{2k-l)\. Let |A| = Y.i^u and length /(A). Thus, /(A) is the 
number of non zero elements of the partition A, i.e. Aj = for j > /(A). Let mj{X) be 
the number of J 's in the partition A, so that |A| = 77ii(A) + 2777-2 (A) + 3m3(A) + .. .. There 

are {ii%){rmix)^rnliX)J monomials in mx{xi, . . . ,Xk) ([S], 7.8). Here {^^^x)^^l^x),.J ^^ 
the multinomial coefficient. Since we are working with symmetric functions, it is enough 
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d a 



to compute (|2.32p . i.e. ()2.34p . for one monomial of m\ [-^, • • • , g^j- Therefore, 



d d \ 



k 



equals 
(2.36) 

l(A);Ui(A),m2(A),...y'5^A, ^K.) \^x^■■■^xJ [dxj- ' ' dxl) l}/^''^> 

evaluated at (xi, . . . , Xk) = (1, . . . , 1). We already have the expression for the effect of 
Vandermonde determinant operators in (j2.25p . Therefore by Lemma [2.3t the expression 
([2:36]) equals 
(2.37) 

(rD( .^^^^^l^^ ) ^^^ E Sgn(/x) det (/('^-l)+20-l))(x.)) • 

V/(A);V"^i(A),m2(A),...;9^A.___^^V, Z.^ IM^ V^ ^)J 

The expression (I2.37P is equal to 

(2.38) (, !",,)( ,J^^\,, ) Vsgn(^)det(/('^-^+20-i)+A.)(i)A 

In each summand of (I2.38p . rearrange the rows so as to reverse the effect of //. We get 

^ ^ \l{X)J\7m{X),m2{X),...J J^^ ^' "^ '> 

Here v is /i^^. The expression (j2.39p is 

^^•^°^ \1{X)) Ui(A), m2(A), . . . j J^^ ^'* ((2fc_i_(,_i)_2(j-l)-A.J!J,. " 

Each determinant inside the sum is of the form 

1 



(2-41) det , , , , , , ^ ,x, 

and ^ dj = |A|. In Proposition 12. 5|. we determine a necessary condition for the determi- 
nant (j2.4ip to be non zero. This condition will imply that a large portion of terms in 
(j2.40p are zero. 

Proposition 2.5. Consider the determinant 

1 



(2.42) det, , , , , ^ .^, 

\{2k-l-{i-l)-2{j-l)-di)\,^^ 

Assume that J2i=idi = \X\, with di G 'L-yQ.The determinant (j2.42p is zero if any of 
di, . . . , (ifc-IAI ^'^ '^'^'^ zero. 
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Proof. Let uhe a number between 1 and k such that du is non zero. The u^^ row in the 
matrix is 

1 



d„V. 



^^■^^^ \i2k-l-{u-l)-2{j-l) ^ur-yi<,<k 

Now look at the row which is du rows below the row u in the matrix ()2.42p . Let this be 



row V where v 
(2.44) 



u + du- Row V, 



1 



l)-d„)\ 



is identical to row u if dv is zero. We have a necessary condition for the matrix to have 
a non zero determinant; for every u such that du ^ 0, either du+du is also non zero or 
u + du > k. We look at this cascading process, and see that if we start at a row above the 
row k — |A|, that is if (i„ 7^ for some u < k — |A|, then we cannot go down beyond row 
k since all di add to |A|. Hence we will have two identical rows. We can then conclude 
that we obtain non zero determinants in (I2.42P only when du = for l<'u<fc— |A|. D 

The above proposition tells us that, in (j2.46p all the action takes place in the last |A| 
rows or lower. Thus, let u = {uk, . . . , ui) be a permutation of (Ai, . . . , Afc). Notice that 
we have reversed the order of the subscripts on the u's, starting at k and ending at 1. 
Applying the above proposition, we shall assume Ui = for i > |A|. Note, however, that 
some of the Uj's, with i < |A| can also equal 0. For a given permutation u, let i(u) be 
the smallest positive integer such that Uj = for all i > i{n). Thus, i(u) < |A|. 

Next, any two permutations that have identical non-zero Uj's, i.e. that move around 
the O's, produce the same determinant. For any given way of selecting where the non- 
zero Aj's go, there are {k — ^(A))! ways to move around the remaining zero- valued Aj's. 
Furthermore, permuting identical non-zero Aj's also produces the same determinant. 
For a given permutation, there are mi(A)!m2(A)! . . . ways to move around the identical 
non-zero Aj's. Using the fact that 

k\f /(A) 

l{X)J\mi{X),m2{X), 

and taking into accou 
as 



(2.45) 



(fc-/(A))!mi(A)!m2(A)! 



A:!, 



(2.46) 



kij:' 



he above two 

1 


pa 


ragraphs, expression 
1 


(|2.4()|) 

1 


(2fc-l)! 


(2fe-3)! 


. . . -^, 
1 


(2fc-2)! 

1 


(2fe-4)! 
1 


••• 0! 


(fe+*(u))! 


{k-i{u)-2y. 




1 




1 




(fc+i(u)-l-ni(„ 
1 


))' 


(fc+J(u)-3-«,(„))! 
1 


... 


(fc-ui)! 


(A:-ui-2)! 
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There are k — i{u) rows above the horizontal dashed hne and i{u) rows below the dotted 
line. The sum is over distinct permutations {uk, ■ ■ ■ , ui) of (Ai, . . . , Afc), satisfying m = 
for i > |A|. Note that, in order for a given permutation u to appear in the sum, we 
require that k > i(u). 

We may also reduce the number of terms in the sum by excluding matrices where two 
or more rows of the matrix are identical. The ' on the sum indicates that such terms 
have been excluded from the sum. 

Now consider one specific term in the sum (j2.46p . As in the calculation of the leading 
coefficient, multiply its i row by {2k — i)l and its j column by l/(2(j — 1))!. This 
enables us to write the determinant in a term of (I2.46P as a product of a known quantity 
and a determinant of binomial coefficients, 



(2.47) 



n,ti(2(j-l))! 



/2k-l\ 

V > 
/2k-2\ 

V ) 






^k+i{u)\ 
/fe+i(u)-l-u,(u)-j 



/fc+i(u)\ 

/fc+i(u)-l-Mi(u)\ 



f2k-l\ 
\2k-2) 
j2k-2\ 
\2k-2) 



(k+i(u)\ 
\ 2k-2 



( 



fc+^{u)-l-u,(^)^ 
2k-2 I 



t^') 



(k-ui\ 
\2k-2) 



X {k + i(u) 



CT) ( 

Here (x)„ is the falling factorial x{x — 1) . . . (x — n + 1). The last factor, the product 
of falling factorials, is a polynomial of degree |A| in k. Expressions (I2.46P and (I2.47P 
for the lower terms are the analogue of (12.28P for the leading coefficient. The difference 
is the presence of ui, . . . ,Uj(u) in the determinant and the appearance of the product 
of falling factorials. The latter are accounted for by the fact that the {2k — i)! is not 
entirely cancelled by the numerator of the binomial coefficients in the last i(u) rows. 

We study the above determinant in the next section. To apply the formulas of that 
section, we require ui > U2 > ■ ■ ., which does not typically hold for the terms in the sum 
of (j2.46p . However, by swapping adjacent rows, we can arrange that these inequalities 
hold. More precisely, say that Um < Um+i- We can assume that, in fact, Um + 2 < Um+i 
since if Um + 1 = Um+i then the ?n,-th and m + 1-st rows from the bottom of the matrix 
coincide, and such terms are excluded from (j2.46p since the determinant in such cases 
is 0. 

Consider what happens when we swap the ?Ti--th and m + 1-st rows from the bottom. 
The binomial coefficient (''+™~j^^''™) gets switched with {^^"^~^^+^) at a cost of a sign 
change to the determinant. On the other hand, the new determinant is of the same 
form, but with u replaced by u', where u'j = Uj for all j, except for n^ = Um+i — 1 and 



u'^+i = Um + 1. Thus we have reversed the inequality, i.e. u'^ > 



u'm+i- Notice also that 



this swapping also satisfies Y^u'j = Y^Uj = |A|. 

Therefore, continuing in this fashion, any given determinant in the sum in (I2.47P is 
equal, up to a power of —1, to the same kind of determinant but with u replaced by, 
say, a(u), where a is a partition of |A|, i.e. with qi > Q2 > . . . > 0. Let the power of 



20 GOULDEN, HUYNH, RISHIKESH, AND RUBINSTEIN 

— 1 introduced by the row swaps that take u to a(u) be denoted by n(a). Thus, a given 
determinant in (j2.47p is equal, on performing the row swaps, to 



(2.48) (-l)"("^+(')l?«(u)(A:), 



k\ 



where D is the determinant defined in ()1.23p . The extra (— 1)^2J arises because the 

columns of D in (jl.23p are in the reverse order from the determinants in (j2.47p . 

Therefore, returning to equations (j2.6p . (j2.34p . and (j2.32p . we have, on simplifying, 

fc(fc+i) 
that the coefficient c±{r,k) of x 2 ''in Q±{k,x) can be expressed as: 

X E ^A (^) E '(-l)"^"^^a(u)(A;) X (A; + i(u) - 1)„,(^, {k + i(u) - 2),^(^,_^ • • • {kU. 

\X\=r u 

In Theorem 11.21 and Corollary 11.31 we show that, for k > max(/(a), ai), 

(2.50) D^{k) = 2(2)-P^{k), 

where Pa{k) is a polynomial in k of degree |q|. Theorem 11.21 also gives a formula for 
determining the polynomials P. Hence 

(2-51) c^ir,k) = ("l n (ifTyr) i:fx(k)N,ik), 

where 

(2.52) N^{k) = E '(-l)"^"^^a(u)(fc) X (A; + i{u) - l)^^^^, {k + .(u) - 2)„^(^,_^ • • • [k^. 

u 

The sum is over distinct permutations u = (uk, . . . ,ui) formed from the partition A by 
appending O's if necessary. Furthermore, we have restricted to permutations such that 
^(u) < |A|, and have also excluded u where the corresponding matrix has any identical 
rows. Finally, for a given u to appear in the sum we have assumed that i(u) < k, and 
to apply (|2.50p . we also required that k > max(/(a(u)), a(u)i). 

We show that the latter assumption can be removed. First note that i(u) = /(a(u)), 
because our swapping procedure that replaces a given u with u' has i{u') = i{u). Next, 
Corollary 11.41 tells us that Paik) vanishes for ai < k < l{a) — 1. Furthermore, 

(2.53) {k + i(u) - 1)„^(^, {k + i(u) - 2)„^(„j_^ • • • {kU 

vanishes if < A; < ai as can be seen by examining the factor associated to ai: let 
Uj be the term that, under our swapping procedure, gets swapped down to ai. The 
corresponding falling factorial is {k + j — l)uy But ai = Uj — {j — 1), because Uj gets 
moved down j — 1 rows to the bottom row. Therefore, 

(2.54) {k + j- 1)„^, = {k + j- l)a,+j-i = ik + j-l){k+j-2)...{k-ai + l) 
which is divisible by 

(2.55) k{k-l)...{k-ai + l). 
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Thus, we have shown that 



(2.56) 



Paiu) (k) x{k + i{u)- 1)„^(^, {k + i(u) - 2)u^^^^_, • • • {k)u. 



vanishes for < A; < max(/(Q(u)), a(u)i). We can, therefore, ignore, in (j2.52p . the condi- 
tion that k > max(/(Q(u)), a(u)i), since, in including terms with /c < max(/(a(u)), a(u)i), 
the corresponding summand in ()2.52p vanishes. 

Hence, Nx{k) is given by a sum over a fixed, i.e. depending only on A but not on 
k, number of terms u. Each term is a polynomial of degree 2|A| in k, thus Nx{k) is a 
polynomial in k of degree < 2|A|. 

This completes the proof of Theorem 11.11 

As an example. We compute A'"(2,i,i)(A;) using ()2.52p . In this case, the sum (|2.52p is 
over the 12 distinct permutations of (2,1,1,0). We can truncate at 4 terms because 
|A| =4, andi(u) < |A|. Of these 12 permutations, only (2,1,1,0), (0,2,1,1), (1,0,2,1), 
and (1, 1,0,2) give non zero determinants. The sign (— 1)"(") is 1 for (2, 1, 1,0), and -1 
for the rest. We have 



Ni2,i,i){k) = P(2,i,i) {k)2{k + l)i(A: + 2)i - P(i,i,i,i) {k + l)2{k + 2)i(A; + 3)i 
- ^(1,1,1,1) (^)i(^ + 2)2(A: + 3)i - P(i,i,i,i) {k)i{k + l)i{k + 3)2 



1 



1 



{k - 3)ik - 2){k - l){k + i) 



-{k -2){k + 3){k'^ + k-4)xk{k-l){k + l){k + 2) ^^ , 
o 24 

X {{k + l)k{k + 2){k + 3) + k{k + 2){k + l){k + 3) + k{k + l){k + 3){k + 2)) 

(2.57) 

=k{k - l){k - 2){k + 3){k + 2){k + 1) 



Having shown that N\{k) is a polynomial in k of degree < 2|A|, we can determine it 
either using formula (12.52P and the formulas in Theorem 11.21 for the polynomials P, or 
else by evaluating (J2.46P for 2|A| + 1 values of k and applying polynomial interpolation. 
More specifically, we can work back from (j2.46p to ()2.32p . and divide by ||- ni=o (k+j)\ 
to get the formula 
(2.58) 



Nx{k) 



^\'^'-'^'\\x\'-^(k + j)\^, 



Y) 



j=0 



(2j)! 



(2fc-l)! 
(2fc-2)! 



(fc+j(u))! 



(2fc-3)! 
(2fc-4)! 



(fc-J(u)-2)! 



1 1 

{k-uiy. (A;-ui-2)! 



1! 
J_ 
0! 



1 1 

(A:+i(u)-l-ni(u))! (fc+i(u)-3-Ui(u))! 
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This formula can be used for a specific choice of A and several values of k to create a 
table of values of Nx{k) to which polynomial interpolation can be applied. Table [T] lists 
the polynomials Nx{k) for all |A| < 7. 

3. The coefficients 6^(A;) 

In order to compute the multivariate Taylor expansion of (jS.lSp . we consider the 
series expansion of its logarithm. We first examine the arithmetic product, and let 

oo 

(3.1) log{Ak{zi, ..., Zk)/ak) =: XI XI Bx{k)mx{z). 

r=l \X\=r 

We start the sum at r = 1 because the division by a^ makes the constant term 0. Now, 
the Ihs is symmetric in the Zj's, and we can find Bx{k) by applying 

(3.2) 



Ai!A2!...A,!az^5z2' '" d 



Z, 



where I = /(A), and setting zi = . . . = z^ = 0. Since the partial derivatives do not 
involve z^+i, . . . ,Zk we can set these to before the differentiation. Thus, by (jl.lOp . 
Bx{k) is equal to ()3.2p applied to 

(3.3) -log(a,) + X E l°g(l-3T^)+ E(^-01og(l ^ 






.4{h{^-^rHr\u 



1 \ V. 1 \' '^^ \ 1 \ , / 1 



1 + ^— 1 + — +- -log(l + - 

p2+'0 \ P^J I Pj ^ P 



evaluated at zi = ... = zi = 0. Likewise, we can find the coefficients of the expansions 

1 fc oo 

(3.4) --Xlog^(l/2 + ^„a)=:X E ftik)^xi^), 

j=l r=l\X\=r 

where a = 1 for /~ and for /"*", and of 

oo 

(3.5) X log (C(l + Zi + Zj){zi + Zj)) =: E E 9x{k)mx{z), 

l<i<j<k r=l \X\=r 

by applying ()3.2p . at zi = . . . = z; = 0, to 

1 ' 

(3.6) --J2logX{l/2 + z„a), 

and to 

(3.7) X l°g (C(l + ^i + ^d^^i + ^i)) + E (^ - log (C(l + Zi)zi) , 

respectively. 



LOWER ORDER TERMS 23 



Next, by composing the three series expansions ()3.ip . (I3.4p . ()3.5p with the series for 
the exponential function, we can derive formulas for the coefficients b^{k). Example 
formulas, for b,-^Jk) and b,-^ i){k), are displayed in the introduction. 

To obtain numerical approximations to b^ (k) for specific choices of k and A one 
needs to compute infinite sums over primes where the summand is a rational function 
of p^'^ times log(p)' '. This can be achieved to high precision using Mobius inversion as 
described, in the context of the moment polynomials of the Riemann zeta function, in 
Section 4.1 of [CFKRS2] . In this fashion, and using (ll.lSp . we computed the values of 
c±{r, k), for r < 10 and k < 9, given in Tables [3] and IH 

4. Determinant of a matrix of binomial coefficients 
Proof of Theorem \1.2[ We shall first prove ()1.27p , and use it to prove (|1.26p . 

Proof of (fOTl) : 

For a fc-tuple (ai, . . . , Ofc) and x = (xi, . . . , Xfc), let x" denote the monomial x"^ . . . x'^'' . 
For a partition A of length less than or equal to k, x^ can be defined by appending zeros 
after the positive elements of A to make it a /c-tuple. 

Reversing the rows of the matrix in (jl.23p , we see that 



(") °^w -(-')<=' ^-(C '2;: 2;'-)) 



l<i,j<fc 



The (z, j)th entry of the matrix in (14. ip can be written using the coefficient operator 
defined in (ll.25p . Let x = {xi, . . . , x^). Then 

(4.2) D,{k) = {-lpdet{[xf-'^]{l + x,)'+'-'-'^)^^^^^^^. 

Noticing that column j only involves Xj, we can move all the [x- ■'] in front of the 
determinant to get 



{-l)(')[x^^'^]det ({1 + Xj)''+'-^-^^) 



l<i,j<k 



(4.3) = (-l)b)[x2^^-] det ((1 + Xj)-('=-^+^')) n(l + xif'''\ 

1=1 

The determinant in (j4.2p can be written in terms of a^^, and Sx defined in p.53p and 

dLlSD, 

(4.4) D,{k) = (-1)© [x2^^]a,+,,(Tqk' • • • ' TT^) 11(1 + xif'-' 

1=1 

k 

= i-lp [x''^]asdT^^, • • • > TT^)^a(t^, • • • > IT^) 11(1 + xir~\ 

1=1 
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But ()1.56p gives as^{xi, . . . ,Xk) explicitly. Hence 



,1 1 ^_ TT / 1 1 \_ TT C^ + Xj) - {1 + Xi) 



OfeVl+^l,---,l+^,; ,li^,Al + ^i 1 + xJ ,J}^, {l + Xj){l+X^) 



(4.5) 



Ui<i<j<kixj - Xi) _ (-1) (2) 05^,(2;) 



n.tiii + x,)'^-! n'=i(i + ^,)'=-^" 

Using (j4.5p in ()4.4p . we have 

k 

(4.6) Z),(A;) = [x2^^]a4(xi,...,Xfe)5A(TT^,...,Tqi^) 11(1 + ^0'- 

We shall now express sx as a coefficient in a polynomial which is easier to work with. 
The dual Jacobi-Trudi identity, (jl.59p . gives 



(4.7) sx = det (e^^-.+j-; 



l<i,jf<n ' 



where {fii, . . . , /i„) is the conjugate partition of A, and n = /(//). 
Expanding the determinant, we get 

n 

(4.8) sx= ^ sgn{a)Y[e^^^,+^(^i-j{x). 

From (I1.46J) . we have Cr = [t'^]E{t). We rewrite ()4.8p using this notation. Let t 
(ti,...,t„). Then 

n 

sx{x)= Y. sgn(a)[t5'^-^+'^(^)...C-"+'^W]ni?(i.) 

(7£Sn 4 = 1 

n / " \ 

Next, pull out nt^ from the sum, and multiply and divide by JJtf, to get 

n n 

i=l l<i<i<n i=l 

n 



(4.9) =r+'-]asMl[Eiti)- 

i=l 

Here we have also used the Vandermonde determinant (up to (—1)^2''): 

n 

(4.10) E^g-(-)nc^^'^= n (t^-t^y 

ct£S„ «=1 l<i<j<n 
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We have expressed sx{x) as a coefficient in a polynomial. Substituting ()4.9p for sx in 
()4.6p . and using the product form of E{t), (jl.46p . we have 



^A(A;) = [i^+'"x2^^]a4(x)a5„(t) 



n / fc 



n na+T^ 



.4 = 1 \l = l 

k / 



xi ' 



Wil + xi) 



[t^^+'-x^''']as,{x)as^{t) J] (1 + ^O'"" 11(1 + ^' + t,) 



A: 



n 



(4.11) 



[i^+^"x2^^]a,,(x)a,„(t) n ( (1 + ^0'-" 11(1 + if^) I 11(1 + t,)^ 



1=1 



1=1 



Applying the dual Cauchy identity ()1.6ip to the double product on the rhs above gives 



(4.12) n (l + ^0'-"n(l + Tf^) =E^a(xi,...,x,)sv(1,...,1,t^,...,T^). 
1=1 \ i=l I \ 

The number of I's in the second factor on the right hand side of (j4.12p is A: — n. Recall 
from (jl.55p that a^sx = as+x- Hence (|4.1ip equals 

n 

(4.13) [tf^+^'^][x''^>^]as^{t)Y[il+U)'Y.^MS,{xi,...,Xk)syil,...,l,j^^,...,j^). 

i=l A 

The monomial x >" occurs in the sum in (j4.13p only when A = (5^. The coefficient of 
x^*^* in a2(5fe(xi, • • • ,Xk) is 1. Simplifying (|4.13p . we have 



(4.14) Dx{k) = r^'-]as„it) n(l + U)'ss,il, . . . , 1, i^: 



l + tn 



i=l 



Note that we have used 6'f^ = 6k- Applying the formula ss^. in (jl.57p . we have 
(4.15) 



^4(i,---,i,TTir'---'TTt;:) = 2^^ ^11(1 + 1 

i=l 



+ L 



k—n 



l<S<nVl + i.^l 



+ t, 



The 2v 2 J comes from pairing, in applying (jl.57p . the k — n Vs. The middle factor 
arises from matching each 1/(1 + tj) with k — n I's, and the last factor from matching 
all pairs of distinct 1/(1 + tj), 1/(1 + tj). Substituting (j4.15p into (|4.14p . and collecting 
the powers of (1 + tj) gives 



(4.16) Dxik) = r+'"]as,,{t) 2(V) m(i + t.)(2 + t, 



,k—ri 



\i=l 



n (2 + ti + tj). 
l<i<j<n 
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Substituting Zi = tj/2, and collecting powers of 2 (note that ( ^^) + {k — n)n+[^) = (2)), 
we get 



■n 

(4.17) Dx{k) = [z^+^"]as„{z)2(')H2)~\^^+^-\ [](l + 2z,)(l + 



z,:)'-" 



X n (i+zi+zj). 

l<i<j<n 

Here we have also used as^{t) = as^{z)2\2) , Since |(5„| = (2), and |^| = |A|, this proves 
(fOTD . 

Proof of (fT^ : 

We now use (|1.27p to prove (|1.26|) . As above, let z = (zi, . . . , Zn). 

Since Schur symmetric functions form a Z-basis for the ring of symmetric functions, 
the coefficient of s^ of a symmetric function F is well defined. We denote this coefficient 
by [s^]F. 

For a symmetric polynomial F(z) in n-variables, and a partition 7 with length at 
most n, we have 

(4.18) [s,iz)]Fiz) = [z'y+'"]as„iz)F{z). 
This can be seen by writing F{z) in terms of our Schur basis 

(4.19) F{z)=J2v^s^{z). 

7 
We wish to find the coefficient v^. Multiplying by as„{z) and using (I1.55P gives 

(4.20) as^{z)F{z) = ^-^^.a^+^Jz). 

7 

Now, the monomials in 0^+5^(2;) are all distinct, and distinct from the monomials in 
0"y'+5„ for any different partition 7' of length at most n. Furthermore, z'^'^ " appears in 
a-y+s„{z) with coefficient 1, coming from the main diagonal of (ll.53p . Thus, v^ is equal 
to the coefficient of z'^'^ "- in as„{z)F[z). 
Therefore we can rewrite ()1.27p as 

(4.21) D,ik) = 2('^y^%,{z)] n fTT^lTrTSlna + 2-0(1 + ^^)'"'- 

We shall work with the ring of symmetric functions A instead of the ring of symmetric 
polynomials in n variables A„. The right hand side of (j4.2ip equals 



,4.22) 2(^)H.,„(.„ n^ {^,^±1^) n(l . 2.,(: . .)-. 

Note that in (j4.22p . we are looking at elements in the ring of symmetric functions, A, 
i.e. as a product involving a countable number of variables zi,Z2, ■ ■ ., whereas in (14.2ip . 
we were considering the elements in the ring of symmetric polynomials in n variables, 
A„.. 
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Applying w, and using ()1.64p we obtain 
(4.23) 



We use the fact that exp(log(l + «)) = 1 + m to write the argument of u as formal power 
series: 



i>l i>l 



(4.24) = 2(^)-l^l [.,(,)] a;f exp^ i (- E ( .I ( ") E ^V 



a+b„fl+c 
J 



-(-2)"Ey?-(^-i)(-irEy? 

i>l j>l 

We can rewrite the argument of uj in (I4.24p using power symmetric functions; 

(4.25) i?A(A;) = 2©-l^l[s,(y)]a;(expE-(- E 

(T) (T) ^^P-+^P-+- - P^-+b+c) - {-'^TPa -{k- l){-lTPa\ j 

In (|1.64p we have seen that uj{pa) = {—lY~^Pa- This gives 

(4.26) D,{k) = 2(^)-l^l [s,{y)] exp E ' f " E f T) {"") 

a>0 " V 6,c>0 V " / V '^ / 
{-lf''^^^'Upa+hPa+c+P2a+h+c) + 2>a + (A: " l)pa 



,4.7, .,©-,.,,„„ n(i-^r^Jf^)na-%.)-(>-.)- 



If we isolate factors corresponding to i = j in the first product, we are able to cancel 
some factors in the second product. Simplifying, we get 



..8, «.,.,. 2ffl-N,,(,„n(.-j,3^fh;;)) no-..)- 



fc-i 
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To calculate the coefficient of sx in (I4.28p . we only have to look at the projection in 
any A^ such that m is greater than or equal to /(A). We can choose it to be equal to 
1{X) (also equals /ii). Let m = /(A). Then 



(\ m 



l<i<j<m \ ^^ yi)(l -%■)/ j=i 

which is equal to 



(4.30) 2('^y\%xiy)] n ii-m-yj)Ei^-y^r'~"'- 

l<i<j<m i=l 

Another application of (I4.18P proves (ll.26p . 



D 



Proof of Corollary \1.3[ It is immediate from ()1.26p or ()1.27p that P\{k) is a polynomial 
in k with integer values at integers of degree at most |A|. We will show that it is in fact 
of degree |A| and determine its leading coefficient. 

From ()1.26p . the highest power of k occurs when we pick as many powers of yi as 
possible from the last product. This happens when none of the yi are picked from 
{1 — yi — yj). Note that 

(4.31) (1-2/)-'=-™ = l+(fc+m)2/+ ^^ + "^^^^ + "^ + ^^ 2/^+. . . = f:ik+m+j-l), y^/j\ 

The coefficient of the highest power of k that appears in the j'-th term of this Taylor 
series is 1/j!. Thus, the coefficient of /c'^' in Px(/c) equals 

[2/1^^+™-^ . . . y^™] [ n (y^ - y}) ) exp(yi + • • • + 2/m) 

creSm \i=l J 

= E sgn(f7)— — 7T77T =det(l/(Ai-i+j)!)mxm 

J^^ {\^-^ + <J{^))\ 



(4.32) 



. ^l<^<i<m(A^-AJ■-^+j) ^ X^(l) 

Y{i<i<m{h + m - i)\ |A|! 



where x^{^) is the degree of the irreducible representation of S\\\ indexed by A. See 
example 6 in Chapter 1.7 of [M] for the last two equalities. 

D 

Proof of Corollary \1.4\ We use equation (|1.27p which gives a formula for Px{k). As part 
of the process of identifying the coefficient of zf^ "- ... z^" in that formula, we focus on 
the coefficient of Zj*^ "" . Now, /xi = 1{X), and n = Ai, hence iJi+n — l = 1{X) + Ai — 1. 
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When we expand ()1.27p . some of the powers of zi come from the factor (1 + zi)*^^^^, 
and the rest from 

(4.33) n izi-z,)il + zi+Zj){l + 2zi). 

l<j<Xi 

Consider the terms arising from taking a z^ from the above. Notice that ()4.33p is 
a polynomial in zi of degree 2Ai — 1, and thus < j < 2Ai — 1. The remaining 
1{X) + Ai — 1 — j powers of zi come from expanding (1 + zi)^~''^^ using the binomial 
theorem, so that the term associated with a particular choice of j is divisible by 

,, ,. ( A: - Ai \ (fc - Ai)(fc - Ai - 1) ■ ■ ■ (fc - 2Ai - /(A) + 2 + j) 

^ ■ ' yiiX) + Xi-l-jJ (/(A) + Ai-l-i)! 

For all < j < 2Ai — 1, this is divisible by 

(4.35) {k - Xi){k -Xi-l)...{k- 1{X) + 1). 

The coefficient of z^ = z'^^ " in the expression in (jl.27p is therefore divisible 

by dSSD. Thus, so is the coefficient of zf'+""^ ■■■zH". 

The same analysis applied to ()1.26p . and using (j4.3ip gives that Px{k) is divisible, for 
^A) < Ai, by 

(4.36) {k + /(A)) . . . (A; + Ai - l){k + Ai). 

D 

5. Family of quadratic twists of elliptic curve L-functions 

Here we modify our techniques to the family of L-functions associated to the quadratic 
twists of an elliptic curve over Q. To keep things explicit, we focus on the elliptic curve 
of conductor 11: 

(5.1) Eua:y^ + y = x^-x. 

The L- function of Ena is given by an Euler product of the form 

(5.2) Luis) = ^ _ ^-^_,_i/2 n^ ^ _ a{p)p-s-i/2 +p-2.' 
which can be expanded into the Dirichlet series 

n=l 

The Dirichlet series above is absolutely convergent in 3f?s > 1. The coefficients a{n) can 
be obtained from the Fourier expansion of the cusp form of weight two and level 11 
given by 

OO OO 

(5.4) Y. «(^)9" = 9 n (1 - 9")'(1 - </"")', 

n=l n=l 

or, alternatively, by counting points on Eua over the fields F„, p prime. 
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The function Lii(s) has analytic continuation to all of C and satisfies the functional 
equation 

(5.5) Luis)=X{s)Ln{l-s), 

where 

r(3/2 - s) / 27r ^ 2s-i 



The L-function associated to a quadratic twist of Ena has a Dirichlet series of the 
form 

(5-7) Ln{s, Xd) = Yl i/2+s ^d{n), 

n=l 

where d is a fundamental discriminant which we further assume satisfies (d, 11) = 1. 
Lii{s,Xd) satisfies the functional equation 

(5.8) Ln{s,Xd) = Xd{-n)\d\^-^'X{s)Ln{l - s,Xd)- 

When considering the moments of Lii{l/2,Xd) we should restrict L{s,Xd) to have an 
even functional equation, i.e. Xrf(~ll) = 1) otherwise L{l/2,Xd) is trivially equal to 
0. In [CFKRSj . d was also restricted to being negative since it allowed them to exploit 
a theorem of Kohnen and Zagier |KZ] to easily gather numerical data for Lii(l/2,Xd) 
with which to check their conjecture. 

When d < 0, Xd(— 1) = — 1, hence, in order to have an even functional equation, we 
require Xd(ll) = ~1) i-e. d = 2,6, 7, 8, 10 mod 11. CFKRS conjectured, see Section 5.3 
of |CFKRS] ■ the asymptotic expansion: 

1 ^ 1 f^ 

(5.9) E L^i{l/2,xdf ^ -r^^ Tk{logt)dt. 

d=2,6,7,8,10 mod 11 

The extra factor of 5/11 on the rhs, compared to (II. 5p . reflects the fact that the sum on 
the left is over 5 out of 11 possible residue classes mod 11. Here, T^ is the polynomial 
of degree k(k — l)/2 given by the fc-fold residue 
(5.10) 

i=i 
where 

k 

(5.11) Rii{zi,...,Zk) = Ak{zi,...,Zk)X{X{l/2 + z,)-^/^ [] C{l + z, + z,), 

i=l i<i<j<k 

and, overloading the notation of Section [TTTl Ak is the Euler product which is absolutely 
convergent for J2^=i \^j\ < I' 

(5.12) Akizi,...,Zk) = l[Ru,p{zi,...,Zk) n (l-rT+IITi- 

P l<i<j<fe 
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with, for p ^ 11, 
(5.13) 



Rii,p — ( 1 H — 
\ p 

and 



l\-^ f I I I ^ 1 ^ 1 

pJ \ p 2 \ -IJ: 1 — a{p)p^^^^j +p^i~22j 11 ]^ -\- a{p)p^^^^3 -\- p" 




k ^ 



(5-14) ^n,n = n3^:^_,_.. 

Note that, although here we are working with the specific elhptic curve Ena, CFKRS' 
recipe provides a similar conjecture for the quadratic twists of any elliptic curve over Q. 
For many examples, see the paper |CPRWj . The only difference is in the conductor, in 
the local factors of A^ for the primes dividing the conductor, and in the allowed residue 
classes (and modulus) for d. 

Next, Tk{x) is a polynomial of degree of k{k — l)/2 given by the A;-fold residue 
(I5.10p . The degree works out smaller compared to Q±{k, x) because the product of zetas 
in (jS.lip involves fewer zetas, i.e. the product over i < j has (2) factors. Therefore, we 
can write 

(5.15) Tfc(x) = co{k)x'^'-'y' + ci(fc)x'^(^-i)/2-i + . . . + c,t^k-i)/2ik). 

Also note that the exponential in (j5.10p has an x rather than x/2. This will impact the 
powers of 2 that enter into our formulas for the coefficients Cr{k). 

To address the poles coming from the zeta-product nC(l + -^j + Zj) we absorb some 
of the factors of A{zi, . . . , z^) = Y\i<i<j<k{zj — Zi){zj + Zj). Thus, 

T.(x) = ^^-^^ j^—,j...JA,{z,,...,z,)X{X{l/2 + z,)-'l^ 

(5.16) 

X n C(l + z^ + z,){z, + z,) ^^''' ■ -^f ^^^11' ■■■'"'^ exp(x J2 ^^)dzi ■ ■ ■ dz,. 

l<i<j<k llj=l ^j j=l 

We overload notation again and set 
(5.17) ak:=Ak{0,...,0) 

and expand 

(5.18) 
1 fc 00 

—Akizi,...,Zk)l[Xil/2+Zj)-^/^ n Cii+ z^+Zj)izi+Zj)=■.Y, E hik)mxiz), 

"'^ i=l l<i<j<k J=0\X\=j 

where, as before, m\{z) is the monomial symmetric function for the partition A. The Ihs 
above is holomorphic in a neighbourhood of zi = . . . = z^ = 0, because the poles from 
the zeta-product J^ C(l + Zj + Zj) are canceled by the product Y\{zi + Zj). We normalize 
by Ofc so that the first coefficient is 1. 
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So (I5.16P becomes 

(-l)^(^-i)/22fe a, r r^ 
T.(x) = (^-y^f ../ E E/'AlfeWl.) 

(5.19) X ^^"^' • • -f ^^Jj; • • • ' "'^ exp(x E ^.)rf^i • • • ^^fe- 

Comparing to equation (j2.6p we notice three differences: the extra 2 in front of the 
integral, the 2k — 1 powers of each Zj, rather than 2k powers, in the denominator, and 
the X rather than x/2 in the exponential. The first two differences are accounted for 
by the fact that the product over zetas in (|2.3p includes i = j, and this introduces, 
from ()2.2p . an extra 2zj, for each j. Therefore, proceeding as in Section [2.21 we get: 

crik) = ^^ (^-)lf ••/ E hik)m,iz) 

\X\=r 

^^^(^^'•••'^^)^(^---'^^)exp(E^.M.i...d... 



yu. 


^u; 




n,ti 


#- 


Ar 


lalogously to 


(|2.49p. we have 




(5. 


21) Cr{k) = 


fc-i 

2'a, n 

i=o 


(2i)! 

{k + j-iy. 





X E^A(fc)E'(-ir^"^^a(u)(A;)x(fc + z(u)-2)„^(^,(A; + i(u)-3)„^,^,_^...(A;-l),,. 

|A|=r 

where, for a partition a, 

^2A; - i - 1 - Ofc-j+i' 



(5.22) i?,(A:) = det 

V V -J / / i<jj'<fc 

The equation for Cr{k) differs from (j2.49p in the power of 2 that appears, and also some 
of the factorials have an extra — 1 in them. The latter comes from the one missing Zj in 
the denominator of (j5.20p compared to (j2.3p . 

Notice that Ea(k) is very similar to Da{k). The only difference is the extra —1 in 
the binomial coefficient. We can relate the two determinants by taking advantage of 
the entries in the first column of the matrix for Ea{k), which are all except for the 
1,1 entry. Assume, for now, that k > max(/(a) + l,ai) (so, in particular, a^ = 0). 
Expanding along the first column, and then reindexing i,j with i + 1, j + 1: 

(5.23) E^(k) = <^et (("=-'- '-"'-» 

(5.24) = detfP^^-\\-^-"('=-^)-^+^ 

VV 2(fc-l)-2j //i<,,<,_i 

(5.25) = Z)„(A; - 1) = 2('2')-l"l X P„(A; - 1) 
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Also note, while we have assumed that k > l{a), Corollary 11.41 tells us that Paik — 1), 
and hence the rhs of (j5.23p . vanishes for ai + 1 < k < 1(a). Furthermore, by the same 
method as was used around (j2.55p . 

(5.26) [k + i(u) - 2)„^,^, {k + t{u) - 3)„^,^,_, ■■■{k- lU 

is divisible by (A; — 1) ... (fc — ai), and thus vanishes for 1 < k < ai. Therefore, writing 

(5.27) cAk) = 2^+C=^^)-a, [] TTT^ 

jJo (k + j-iy. 

xE^A(A;)E'(-ir^"^^-(u)(A:-l)x(A;+i(u)-2)„^(^,(A;+.(u)-3)„^(^j_^---(/c-l),,, 

|A|=r u 

we can replace the requirement that k > max(/(a) + l,ai) with A; > 0. Finally, us- 
ing (|2.52p . we get, for A; > 0, 

(5.28) crik) = 2'=+('^')-a, J] ..i^^!!.,, E bx{k)N,{k - 1). 

For example, the r = term equals 

A-n ^_^ (2i)! 



(5.29) co{k) = 2''HV)akl[ 



One can verify inductively that this matches the leading term as described in (1.5.26) 
of jCFKRSj : 

(5.30) 2^+('2^) f\ , ^^^'^' „ = 2('=+i)'=/2 TT -A-. 

;=o(^+j-i)! ;=o(2j)! 

One should also pay attention here that the Taylor coefficients b\{k), and also a^, 
depend on the underlying elliptic curve Ena and its a(p)'s. While we can derive similar 
formulas for bx{k) as for quadratic Dirichlet L-functions (see the examples (ll.2ip . ()1.22p ). 
in order to accelerate their numerical evaluation we would need to use the symmetric 
power L-functions associated to the L- function Lii(s). Acceleration for bx{k) will be 
discussed in a forthcoming paper of Alderson and Rubinstein |AR2] . 
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A 


Nxik)/rx{k) 


rx{k) 


1 


k + 1 


{k)i 


1,1 


{k + 2){k + l) 


{k)2/2 


2 





ik)i 


[1,1,1 


{k + 3){k + 2){k + I) 


{kh/Q 


2,1 


{k + 2){k + l) 


{k)2 


3 


-{k-l){k + 2){k + l) 


ik)i 


1,1,1,1 


ik + A){k + 3){k + 2){k + l) 


ik)4/24 


[2,1,1 


2{k + 3){k + 2){k + l) 


ik)3/2 


2,2 





{k)2/2 


3,1 


-(A:-2)(A; + 3)(A: + 2)(A: + 1) 


{k)2 


4 





ik)i 


1,1,1,1,1 


{k + 5)ik + 4){k + 3){k + 2){k + l) 


{k)5/120 


2,1,1,1 


3{k + 4){k + 3){k + 2){k + l) 


{k)4/6 


[2,2,1 


A{k + 3){k + 2){k + l) 


{khl2 


[3, 1, 1 


-{k - 3){k + A){k + 3){k + 2){k + 1) 


{khl2 


3,2 


-2(A;-2)(A: + 3)(A; + 2)(A: + 1) 


{k)2 


4,1 


-2(A;-2)(A: + 3)(A; + 2)(A: + 1) 


{k)2 


5 


2{k - l){k - 2){k + 3){k + 2){k + 1) 


{k)i 


1,1,1,1,1,1 


(k + 6)(A; + 5)(A; + 4)(A; + 3){k + 2){k + 1) 


{k)e/720 


2,1,1,1,1 


4:{k + 5){k + 4:){k + 3){k + 2){k + 1) 


{kh/24 


2,2,1,1 


10(/c + 4)(A; + 3)(A: + 2)(A: + l) 


(fc)4/4 


[2,2,2 





(fc)3/6 


3,1,1,1 


-{k-A){k + 5){k + A){k + 3){k + 2){k + 1) 


ik)4/6 


[3,2,1 


-{k + 3){k + 2){k + l){3k'' + 3k- 40) 


ik)3 


3,3 


{k -2){k- 4)(A; + 5)(A; + 3){k + 2){k + 1) 


ik)2/2 


[4, 1, 1 


-4:{k + 3){k + 2){k + l){k' + k- 10) 


{khl2 


4,2 





{k)2 


5,1 


2{k - 2){k + 3){k + 2){k + l)(/c^ + A; - 10) 


{k)2 


6 





ik)i 


[1,1,1,1,1,1,1 


{k + 7){k + Q){k + 5){k + 4){k + 3){k + 2){k + 1) 


(A:)7/5040 


2,1,1,1,1,1 


f>{k + Q){k + 5){k + 4){k + 3){k + 2){k + 1) 


(fc)6/120 


2,2,1,1,1 


18(A; + 5){k + 4){k + 3){k + 2){k + 1) 


(A;)5/12 


2,2,2,1 


3{){k + A){k + 3){k + 2){k + l) 


{k)4/Q 


3,1,1,1,1 


-{k -5){k + 6){k + 5){k + A){k + 3){k + 2){k + 1) 


ik),/2i 


3,2,1,1 


-2{k + 4)(A: + 3){k + 2){k + 1)(2A:^ + 2k- 45) 


{hU/2 


[3,2,2 


-10(A: -3){k + 4){k + 3){k + 2){k + 1) 


ikh/2 


[3,3,1 


{k -3){k- 5){k + 6){k + 4){k + 3){k + 2){k + 1) 


{khl2 


4,1,1,1 


-6(A; + A){k + 3){k + 2){k + 1)(A:^ + k - Ih) 


{k),/Q 


[4,2,1 


-10(A: -3){k + A){k + 3){k + 2){k + 1) 


«3 


4,3 


b{k -2){k- 3){k + A){k + 3){k + 2){k + 1) 


^2 


[5,1,1 


2{k - 3){k + A){k + 3){k + 2){k + l){k^ + /c - 15) 


{kh/2 


5,2 


b{k -2){k- 3){k + A){k + 3){k + 2){k + 1) 


{k)2 


6,1 


b{k -2){k- 3){k + i){k + 3){k + 2){k + 1) 


{k)2 


7 


-5{k -l){k- 2){k -3){k + A){k + 3){k + 2){k + 1) 


ik)i 



Table 1. We display the polynomials Nx{k), for all |A| < 7. Be- 
cause each monomial of mx^z) contributes the same to (|2.6|) . Nx{k) 

has, as a factor, the polynomial: rx{k) := {il\)){m^(^x)!ml{x),..) = 
(A:)i(A)/(mi(A)!m2(A)!...), where {k)m = k{k -1) . . .{k -m + 1). There- 
fore, rather than display Nx{k), here we list Nx{k)/rx{k). 
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Px{k) 



Partition 



[1 

[2; 

[1,1 

[3; 

[2,1 



[1,1. 

[3 

[2 

[2,1. 

[1,1,1. 

[4 
[3 

[3,1. 

[2,2, 

[2,1,1. 
[1,1,1,1 

[5 
[4 

[4,1. 
[3 

[3,2, 

[3,1,1 

[2,2, 

[2,2,1, 

[2,1,1,1 

[1,1,1,1,1 

[6 
[5 

[5,1. 
[4 

[4,2, 

[4,1,1. 

[3,3, 

[3,2, 

[3,2,1 

[3,1,1,1 

[2,2,2 

[2,2,1,1 

[2,1,1,1,1 

[1,1,1,1,1,1 



k + l 

{l/2){k + l){k + 2) 

{l/2){k-l){k + 2) 

(l/6)(fe + l)(fc + 2)(fc + 3) 

(l/3)(A: + 2)(A;2 + A;-3) 

(l/6)(fe-2)(fe-l)(fe + 3) 

{l/2A){k + l){k + 2){k + ?,){k + A) 

(l/8)(A: + 2)(A: + 3)(fe2 + A:-4) 
(l/12)(fe-2)(A: + l)(fe + 2)(fe + 3) 

(l/8)(A:-2)(fe + 3)(fc2 + A:-4) 

(1/24)(A: - 3)(A; - 2){k - l){k + 4) 

(1/120)(A: + 1)(A: + 2){k + 3)(A: + A){k + 5) 

(1/30)(A: + 2){k + 3){k + 4){k'^ + k- 5) 

(1/24)(A: + 1)(A: + 2){k + 3){k^ + k - 8) 

(l/20)(fe + 3){k^ + 2k^ - llk^ - I2k + 40) 

(l/24)(fe - 2)ik + 1)(A: + 3){k^ + k - 8) 

(1/30)(A: -3)ik- 2){k + A){k'^ + k - 5) 

(1/120)(A; - A){k - 3){k - 2){k - l){k + 5) 

(1/720)(A; + l){k + 2){k + 3){k + A){k + h){k + 6) 

(1/144)(A: -2){k + 2){k + 4){k + 5){k + 3f 

(1/80)(A: + 1)(A; + 2)(A: + 3){k + i){k^ + k - W) 

(1/72)(A; + 3){k + i){k^ + 2k^ - 13k^ - Uk + 60) 

(l/144)(fe - 3)(A; + l)ik + 3){k + i){k + 2)^ 

(1/45)(A; + l){k + 3){k^ + 2^ _ iQp _ nk + 75) 

(1/72) (A; -3){k + 4:){k^ + 2P - 13P - 14fe + 60) 

(1/144) (A; -3){k- 2){k - l){k + 2){k + 3){k + 4) 

(1/80) (A: - 3) (A; - 2) (A; + 1)(A: + 4)(A:2 + A: - 10) 

(1/144)(A: - 4)(A: - 3)(A; + 3)(A: + 5)(A: - 2f 

(1/720)(A; - 5)(A: - 4)(A: - 3)(A; - 2)(A; - 1)(A: + 6) 

(1/5040)(A: + 1)(A: + 2)(A: + 3)(A; + A){k + 5)(A: + 6)(A: + 7) 

(1/840)(A: + 2)(A: + 3)(A; + 4)(A; + 5)(A: + 6)(A:2 + A: - 7) 

(1/360)(A; - 3)(A; + 1)(A: + 2)(A: + 3)(A; + 5)(A; + Af 

(1/336) (A; + 3)(A: + 4) (A: + b){k'^ + 2k^ - I5k^ - IQk + 84) 

(1/360)(A; + 1)(A: + 3)(A: + 4)(fe + 2)^ik^ + k - 15) 

(1/144)(A: + 1)(A: + 3)(A; + 4)(A:^ + 2k^ - Idk"^ - 20k + 108) 

(1/252) (A: + 4:){k^ + 3k^ - 26k'^ - 57k^ + 277A:2 + 306A; - 1260) 

(1/240)(A: - 3)(A; + 1)(A; + 2)(A: + 3)(A: + 4)(A;2 + A; - 10) 

(1/240)(A: - 3)(A; - 1)(A: + 2)(A: + 3)(A: + 4)(A;2 + A: - 10) 

(1/144) (A: - 3)(A: + 1)(A; + 4:){k^ + 2A;3 - 19A;2 - 20A; + 108) 

(1/336) (A; - 4) (A: - 3) (A: + 5)(A;4 + 2k^ - 15k^ - 16k + 84) 

(1/360)(A: - 3)(A; - 2)(A; - 1)(A: + 2)(A: + 4)(A;2 + A; - 15) 

(1/360)(A; - 4)(A: - 2)(A: + 1)(A: + 4)(A: + 5)(A: - 3)^ 

(1/840)(A; - 5)(A: - 4)(A: - 3)(A; - 2)(A: + 6)(A:2 + A: - 7) 

(1/5040)(A: - 6)(A: - 5)(A: - 4)(A; - 3)(A; - 2)(A: - 1)(A: + 7) 



Table 2. Table oi Pxik) 
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C__(r, 1) 



c-(r,2) 



c_(r,3) 




1 
2 
3 

4 
5 
6 



3.522211004995827732e-01 
6.175500336140218316e-01 



1.238375103096108452e-02 
1.807468351186638511e-01 
3.658991414081511628e-01 
-1.398953902867718369e-01 



1.528376099282021425e-05 
8.968276397996084726e-04 
1.701420175947633562e-02 
1.093281830681910732e-01 
1.358556940901993748e-01 
-2.329509111366616925e-01 
4.735303837788046866e-01 



(r,4) 



c-{r, 5) 



c-{r, 6) 



3.158268332443340154e-10 

1 5.062201340608140133e-08 

2 3.252070477914552180e-06 

3 1.065078255299183117e-04 

4 1.865791348720969960e-03 

5 1.658674128885722146e-02 

6 5.985999910494527870e-02 

7 5.231179842747744717e-03 

8 -1.097356193524353096e-01 

9 5.581253300381869842e-01 
10 1.918594095122517496e-01 



6.712517611066278238e-17 
2.341233253582258184e-14 
3.571169234103129887e-12 
3.127118490785452708e-10 
1.734617312939144360e-08 
6.342941105701246722e-07 
1.541064437383931078e-05 
2.441498848686470880e-04 
2.39092828457395691 le-03 
1.275610736275904766e-02 
2.430382016767882944e-02 



1 .0360046454270032766-25 
6.796814066740219201e-23 
2.037808336505920108e-20 
3.698051408075659748e-18 
4.534838798273249707e-16 
3.972866885083416336e-14 
2.563279107875100164e-12 
1.237229229636910631e-10 
4.491515829566301398e-09 
1.222154548508955419e-07 
2.461203700713661380e-06 



(r,7) 



(r,8) 



c_(r,9) 



8.864927187204894781e-37 

1 9.894437508330137269e-34 

2 5.176293026015439716e-31 

3 1.686724585610585967e-28 

4 3.837267516078630273e-26 

5 6.474635477336820480e-24 

6 8.402114103039537077e-22 

7 8.581764459399681586e-20 

8 7.002464589632248733e-18 

9 4.607034349981096374e-16 
10 2.455973970379903840e-14 



3.372009502181036150e-50 
5.951191608649093822e-47 
5.002043249634522587e-44 
2.664702289380503418e-41 
1.010164553397544484e-38 
2.900498887294046119e-36 
6.555588245821587108e-34 
1.196609980002393296e-31 
1.795828629692653400e-29 
2.244368542496810519e-27 
2.357312576663548340e-25 



4.7277357965875261 13e-66 

1.248019487993274422e-62 
1.585820955757896443e-59 
1.291823649274241834e-56 
7.58066062423973821 le-54 
3.413900516458523702e-51 
1.227404779731471396e-48 
3.618608212113140382e-46 
8.916974338520402569e-44 
1.862786263819570034e-41 
3.334524507937658586e-39 



Table 3. The coefficients c^{r, k) of Q^{k). 
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c+(r, 1) 



c+(r,2) 



c+(r,3) 




1 
2 
3 
4 
5 
6 



3.522211004995827732e-01 
-4.889851881547797041e-01 



1.238375103096108452e-02 
6.403273133040673915e-02 
-4.030985462971436450e-01 
8.784723252866324383e-01 



1.528376099282021425e-05 
6.087355322740111135e-04 
5.189536257221761054e-03 
-2.070416696161206729e-02 
-4.836560144295628388e-02 
6.305676273169569246e-01 
-1.231149543676485214 



.(r,4) 



c+(r,5) 



.(r,6) 



3.158268332443340154e-10 

1 4.070002081481211197e-08 

2 1.961035634727995841e-06 

3 4.187933734218812260e-05 

4 3.233832982317403053e-04 

5 -7.264209058002128044e-04 

6 -9.741303115420443803e-03 

7 6.254058547607513341e-02 

8 5.338039400180279170e-02 

9 -1.125787514381924481e+00 
10 2.125417457224375362 



6.712517611066278238e-17 
2.024913313371989448e-14 
2.611003455556346309e-12 
1.870888923760240058e-10 
8.086250862410257040e-09 
2.126496335543600159e-07 
3.194157049041922835e-06 
2.120198748289444789e-05 
-3.390055513847315853e-05 
-7.750613901748660065e-04 
3.339978554290242568e-03 



1 .036004645427003276e-25 
6.113326104276961713e-23 
1.632224321325099403e-20 
2.605311255686981285e-18 
2.766415183453526818e-16 
2.056437432501927988e-14 
1.095709499896029594e-12 
4.206172871179562219e-ll 
1.149109718292255815e-09 
2.154509460431619112e-08 
2.543371224701971233e-07 



_(r,7) 



c+(r,8) 



,(r,9) 



8.864927187204894781e-37 

1 9.114637784804059894e-34 

2 4.370089613567423486e-31 

3 1.297363094463138851e-28 

4 2.670392092372496088e-26 

5 4.043466811338890795e-24 

6 4.663148139710778893e-22 

7 4.183154331210266578e-20 

8 2.954857264190019988e-18 

9 1.652770327042906306e-16 
10 7.319238365079051443e-15 



3.372009502181036150e-50 
5.569826318573164385e-47 
4.368642207198861832e-44 
2. 164658555649376388e-41 
7.604817314362535383e-39 
2.015327809331532264e-36 
4.184593239584908611e-34 
6.980465161514108456e-32 
9.516651650236242059e-30 
1.073015400698217206e-27 
1.008662233782716849e-25 



4.727735796587526113e-66 
1.181 182697783246367e-62 
1.417926553457661234e-59 
1.089051480593133551e-56 
6.012641112088390226e-54 
2.541594397695401893e-51 
8.555207141044511720e-49 
2.354807833463352272e-46 
5.400892227120418237e-44 
1.046573394851932219e-41 
1.731269798305270612e-39 



Table 4. The coefficients c-|-(r, /c) of Q-^{r,k). 



